Predicted Constraints on Cosmic String Tension from Planck and Future CMB 

Polarization Measurements 
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We perform a Fisher matrix calculation of the predicted uncertainties on estimates of the cosmic 
string tension G/j, from upcoming observational data (namely, cosmic microwave background power 
spectra from the Planck satellite and an idealized future polarization experiment). We employ 
simulations that are more general than others commonly used in the literature, leaving the mean 
velocity of strings, correlation length of the string network, and "wiggliness" (which parametrizes 
smaller-scale structure along the strings) as free parameters that can be observationally measured. In 
a new code, StringFast, we implement a method for efficient computation of the Ce spectra induced 
by a network of strings, which is fast enough to be used in Markov Chain Monte Carlo analyses 
of future data. Performing a calculation with the string parameters left free results in projected 
constraints on G/j, that are larger than those obtained by fixing their values a priori, typically by a 
factor of ~2-7. We also find that if G/j, is equal to the current observational maximum, Planck will 
be able to make a confident detection of strings. However, if Gfi is two orders of magnitude smaller, 
even a perfect, lensing-free measurement of polarization power spectra will not be able to detect a 
nonzero string tension at better than 2a confidence. 

PACS numbers: 98.80. Cq, 98.80.Jk, 11.27. +d 



I. INTRODUCTION 

In the last few decades, we have seen the arrival of an 
era of "precision cosmology," as observational data have 
become precise enough to place serious constraints on our 
theories for the large-scale behaviour of the observable 
Universe. These data have led to the formulation of the 
standard model of cosmology (SMC; see, e.g., Ref. [1]), 
in which the Universe began in a hot, dense state and has 
been expanding and cooling for billions of years. How- 
ever, the SMC is far from complete: there are mysterious 
substances, dark matter and dark energy, whose presence 
we can infer from data but have yet to fully understand, 
and there are also proposals for amendments to the model 
that incorporate as yet unobserved phenomena. 

One such proposal is the existence of topological de- 
fects, especially cosmic strings (e.g. [2-4]). At one time, 
strings were thought to play an important role in gener- 
ating structure in the Universe by seeding density per- 
turbations at early times, but observations of the cosmic 
microwave background radiation (CMB) by the COBE 
satellite and later experiments contradicted this pic- 
ture [5-7] , instead favouring structure formation through 
quantum fluctuations that were amplified by a period of 
rapid inflation. 

Nevertheless, cosmic strings have remained a subject 
of active study, since they are a generic prediction of 
superstring-inspired "brane inflation" theories [8] and ap- 
pear as a general feature of inflationary scenarios inspired 
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by grand- unification-scale physics [9]. There are several 
proposed observational signatures of strings, including 
wedge-shaped wakes in 21cm redshift surveys [10, 11], 
spatial correlations between anisotropies in 21cm radi- 
ation and the CMB [12], gravitational waves from the 
decay of string loops [13, 14], and gravitational lensing 
effects [15-18]. Strings would also leave imprints on the 
CMB, such as step discontinuities in the temperature 
map [19] and B-modes of polarization [20, 21] (although 
the detection of B-modes may not be a "smoking gun" 
for strings, or for any other theory for that matter — see 
Ref. [22] for a recent discussion). 

Because measurements of temperature anisotropies at 
small angular scales become dominated by the effects of 
hot gas in galaxy clusters and clustered infrared-emitting 
galaxies [23] , more precise measurements of polarization 
power spectra should be able to push cosmic string con- 
straints to well below what is possible with temperature 
information alone. Indeed, several ongoing and upcom- 
ing projects (Planck [24], COrE [25], ACTPol [26], SPT- 
Pol [27], BICEP2/Keck [28, 29], POLARBEAR [30], and 
SPIDER [31], among others) will soon measure CMB po- 
larization to unprecedented accuracy, and so in the near 
future we will be able to obtain strict limits on the cosmic 
string content of the Universe. 

Consequently, this work proposes to use simulations of 
CMB polarization (and temperature, to a lesser extent) 
induced by cosmic strings to forecast the extent to which 
strings can be detected by upcoming measurements — in 
particular, by Planck or an idealized future measurement 
of polarization. Research that utilizes simulations of cos- 
mic strings is ongoing [32-37], and our contribution will 
be twofold. 

First, in the spirit of "precision cosmology," we will em- 
ploy a more general model for strings than is commonly 
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in use — namely, the one-scale model of, e.g., Ref. [38], 
but with the mean velocity, correlation length, and "wig- 
gliness" of the simulated network of strings left as ad- 
justable parameters, in addition to the (dimensionless) 
string tension Gfx. Taking these properties as free, and 
in principle measurable from data, instead of fixed a pri- 
ori leads to a weakening of projected constraints on G/i, 
which we calculate using a Fisher matrix approach. We 
also develop a new code for quickly calculating the CMB 
spectra induced by a network of strings, which can be 
applied to the Markov Chain Monte Carlo (MCMC) [39] 
estimation of the one-scale string model parameters, and 
which will be made publicly available for use in other 
studies. 

Second, we will perform forecasts both for a realistic 
observational scenario (CMB measurement by the Planck 
satellite) and a "best-case scenario" situation in terms of 
detecting strings (noiseless measurement of polarization 
across a wide range of angular scales) , incorporating the 
effects of gravitational lensing on the measurement. The 
latter will tell us how well we can ever do in constrain- 
ing the string tension using CMB data, and also give an 
idea as to the range of string tensions we might be able 
to reliably discriminate from a Universe without strings 
(Gfx = 0). 



II. MODELLING COSMIC STRINGS 

The direct simulation of a network of cosmic strings 
without any simplifying assumptions is essentially impos- 
sible, as any such simulation would have to encompass an 
enormous range of physical scales to track the behaviour 
of strings from the early Universe until the present day. 
The standard approach uses the fact that, under some 
basic assumptions about string formation and decay, a 
network of strings evolves toward a scaling solution in 
which average properties of the network (such as the cor- 
relation length) scale in tandem with the expansion of 
the Universe. It is then possible to numerically solve 
either the Nambu or Abelian-Higgs (AH) equations of 
motion for a string, compute an unequal-time correla- 
tor of the stress-energy tensor of a simulated string net- 
work [40] , and use this along with standard perturbative 
techniques to compute the resulting CMB anisotropics. 
However, Nambu and AH simulations generate networks 
with different correlation lengths and mean velocities (see 
Ref. [34] for details). 

Another strategy is to represent a string network as 
a random collection of unconnected, straight segments, 
obtain the stress-energy tensor from these segments, and 
average the resulting CMB spectrum over many realiza- 
tions of the network. This model has been implemented 
in the public code CMBACT [38, 41], which also utilizes 
the equation of state for "wiggly" strings [42, 43], in- 
corporating the uncertainty about small-scale structure 
along the strings via a "wiggliness" parameter, a. The 
code calculates temperature anisotropies, Cj T , as well as 



E- and B-modcs of polarization, Cf and Cf , and the 
cross power spectrum, Cj E , including the effects from 
scalar, vector, and tensor components of the source func- 
tion. We use CMBACT as the basis for the numerical 
calculations in this work. Note that CMBACT does not 
include a detailed treatment of production of small loops 
associated with string crossings, and also neglects the 
subsequent decay of these loops into gravitational radia- 
tion. 

The parameters £ (comoving correlation length of the 
network), v (rms velocity of strings), and a take different 
values in the matter- and radiation-dominated eras, with 
their specific time evolution determined by CMBACT 
primarily from the effects of string expansion and loop 
production [44-46]. However, it is possible to explore a 
wider range of string properties if these three parameters 
are approximated as constant in time, so that they be- 
come free parameters of the model [20, 34], determined 
only by their initial conditions. Along with the string 
tension G/z (written as G\i / c 2 with dimensions included) , 
we then have a model for cosmic strings with four free 
parameters, that opens the door to such activities as 
MCMC estimation of these parameters or forecasts for 
how accurately they can be measured by future observa- 
tions (this work is concerned with the latter endeavour, 
focusing on Gfj, in particular). We illustrate the effects 
of varying these parameters on CMB spectra in Fig. 1. 

It has already been shown that these parameters can be 
chosen to give spectra that match the output of Nambu 
and AH simulations to a reasonable level of accuracy [34] . 
While it is true that future string calculations may be 
able to include a wider range of detailed physical effects, 
we believe that our four-parameter approach at least has 
the virtue of being more realistic that the one-parameter 
assumption that is commonly adopted. 



III. STRINGFAST: A CODE FOR FAST STRING 
CALCULATIONS 

For a single set of string parameters, using CMBACT 
to compute the resulting CMB spectra can take anywhere 
from several hours to several days on a single modern 
CPU. The required time depends on the number of re- 
alizations of the string network, which is directly related 
to the accuracy of the G^s generated as output, since the 
final CiS are averages over all realizations. Clearly, these 
computations must be optimized significantly before ap- 
plications that require large numbers of Ci calculations, 
such as MCMC studies, can become feasible. 



A. Strategy 

It has been proposed that "morphing" techniques bor- 
rowed from computer graphics could be used for efficient 
calculations of CMB spectra [47]. This morphing ap- 
proach defines certain control points based on special fea- 
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FIG. 1: Comparison of CMB power spectra (sum of scalar, 
vector, and tensor modes) when a single string parameter 
is varied while the others are kept constant. In each panel 
the black (dotted) line corresponds to a model with G/j, = 
6.4 x 10" 7 , a = 1.05, v = 0.40, and £ = 0.35. The leftmost 
frames show £ = 0.15 (green, dashed) and £ = 0.5 (blue, dot- 
dashed); the centre frames show v = 0.3 (green, dashed) and 
v = 0.8 (blue, dot-dashed); and the rightmost frames show 
a — 1.0 (green, dashed) and a = 1.5 (blue, dot-dashed). 



tures of the spectra (such as acoustic peaks and troughs), 
interpolates between pre-calculated curves to determine 
where these points would be for a different set of input 
(e.g. cosmological) parameters, and then smoothly trans- 
forms the two adjacent curves to a new curve that pre- 
serves the key features of the original curves. 

The spectra generated by cosmic strings are domi- 
nated by only a small number of key features, such as 
a low-£ bump from rcionization and a larger peak be- 
tween t ?s 500 and 1000 from recombination (see Fig. 2). 
Therefore, it is unnecessary to implement a full morphing 
scheme — rather, the spectra can be described by a small 
number of fitting functions, or, in cases where there is 
nontrivial variation of a spectrum with respect to string 
network parameters, by cubic splines. There is then a 
straightforward strategy for calculating the spectra for 
arbitrary values of the parameters (similar to what is 
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FIG. 2: Scalar (green, dashed), vector (blue, dotted), and ten- 
sor (brown, dot-dashed) modes of CMB power spectra gener- 
ated by strings with Gfi = 6.4 x 10~ 7 , a = 1.05, v = 0.40, and 
^ = 0.35. The sum of inflationary scalar, vector, and tensor 
modes generated by CAMB, using WMAP7 parameters and 
r = 0.36, are also shown (red, solid). 



done in some inflationary CMB codes, such as DASh [48] 
and Pico [49]): 

1. Pre-calculate the spectra on a grid in parameter 
space, with limits based on reasonable priors for the 
parameters and spacing chosen to give acceptable 
accuracy in the approximate curves generated by 
the method. 

2. Fit each prc-calculatcd spectrum with fitting func- 
tions or cubic splines. 

3. Create multi-dimensional splines, as functions of 
the model parameters, for each fit parameter (or, if 
a curve is fit by splines, for each Ci value used in 
that fit). 

4. For the desired set of model parameters, use these 
splines to calculate the corresponding fit param- 
eters (or spline control points), and use these to 
output the resulting spectra. 



4 



It is possible to encapsulate this process into a dis- 
tributable Fortran module, which we dub StringFast. 
From an end-user perspective, the first two steps will 
have been completed beforehand, and files of fit parame- 
ters and/or spline control points will be distributed along 
with the source code. At runtime, one first calls subrou- 
tines which read in these files and assemble the multi- 
dimensional splines into memory. One can then use the 
main routines of the module, which perform step 4, to cal- 
culate a specific CV value as a function of spectrum type 
(TT, EE, or BB, as well as scalar, vector, or tensor), 
£, string model parameters, and optical depth r of the 
background cosmology (see below). The resulting time 
for StringFast to compute a complete (2 < I < 3000) Ce 
curve is less than a second on a modern CPU. 

In addition to the four parameters of the one-scale 
string model, we must account for how the string spec- 
tra change with variation of the optical depth r, since 
this value has an important effect on the height of the 
reionization bump, and also scales the height of the main 
peak by a factor of ~ e~ 2r (as roughly e~ r of the photons 
are scattered on their way from the last-scattering sur- 
face). For the other cosmological parameters, however, 
we assume that variations of their values consistent with 
current observational constraints only have higher-order 
effects (that can be safely neglected) on the C^s generated 
by strings. The code evaluates the string spectra using 
the best-fit parameter values from the WMAP seven-year 
data analysis [50] (although this could be changed by run- 
ning a new grid of spectra). 



B. Accuracy 

For the spectra pre-computcd on the grid in param- 
eter space, we run CMBACT with 12,000 network re- 
alizations. To determine the precision of these spec- 
tra, we choose several sets of sample parameters, gen- 
erate eight sets of 12,000-realization curves, and cal- 
culate the ideal signal-to- noise (S/N) of the maximum 
difference signal. Specifically, we calculate the maxi- 
mum value of (Cf - Cf) /ACf B over 2 < I < 3000 
and over all pairs (A, B) of our eight test curves, where 
ACf B = y/2/(2t+\){Cf- + Cf) (the "noise" being 
merely the cosmic variance of the difference signal) . This 
measures the uncertainty of a simulated curve at any 
point along the curve, as compared to the noise inherent 
in a measurement of the Ces. 

The maximum S/N we find is 1.40, and in fact most 
types of spectra have maximum S/N < 1. Thus, the 
variation in the pre-computed spectra does not rise sig- 
nificantly above the fundamental uncertainty of power 
spectrum observations. Moreover, since the underlying 
model we employ is itself only an approximation of the 
more intricate physics of cosmic strings (ignoring, for ex- 
ample, the details of the production of loops), we are 
confident that this is an adequate level of accuracy for 
the forecasts we wish to undertake, and also for more 
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TABLE I: Attributes of grid used for accuracy tests in Ta- 
ble II. 



Reduced x*~, comparing 
approximate curve to 









"exact 


" curve 






input parameters 


TT 


BB 


EE 


T 


a 


V 




s 


v 


t 


v 


t 


s 


v 


t 


0.087 


1.26 


0.58 


0.19 


5.5 


64.1 


11.0 


10.7 


8.6 


2.0 


23.6 


6.0 


0.088 


1.05 


0.40 


0.35 


3.4 


1.3 


2.9 


10.1 


21.9 


5.2 


8.2 


17.5 


0.090 


1.20 


0.77 


0.27 


6.9 


1.1 


2.0 


3.7 


14.1 


5.3 


3.6 


8.9 


0.090 


1.30 


0.40 


0.23 


9.0 


26.9 


0.3 


3.3 


3.4 


1.1 


5.6 


3.5 


0.091 


1.08 


0.78 


0.44 


9.7 


8.0 


2.2 


8.8 


5.7 


6.6 


7.7 


3.7 


0.091 


1.15 


0.33 


0.42 


3.5 


1.1 


2.4 


2.0 


3.6 


13.5 


1.9 


3.6 


0.092 


1.20 


0.47 


0.28 


5.8 


4.3 


4.9 


6.5 


2.0 


5.8 


9.5 


3.5 


0.094 


1.43 


0.51 


0.33 


1.3 


1.6 


1.3 


0.9 


1.9 


2.2 


0.9 


0.9 



TABLE II: Reduced \ 2 (i- e - xV^OOO) values for approximate 
spectra generated from the specified input parameters, using 
the grid described in Table I. We report results for scalar, 
vector, and tensor modes, denoted by "s," "v," and "t" re- 
spectively. 



general studies of the average properties of a network of 
strings. 

To assess the accuracy of the approximate (off-grid) 
spectra generated by StringFast, we compute a \ 2 statis- 
tic comparing the approximate curves to "exact" (96,000- 
realization) curves generated by CMBACT with the same 
parameter values: 



3000 ( n a 

x 2 = E — 
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£<CMBACTy 
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where the ac e are computed from eight sets of test curves 
and their average. For a selection of test parameters, 
approximate curves are generated using a pre-computcd 
grid described in Tables I, with resulting reduced x 2 s 
given in Table II. Since the C^s depend on G[i through 
the simple analytical relationship Cg oc (G/i) 2 , this pa- 
rameter is not included in the grid. The range of param- 
eters in Table I is sufficient for Fisher matrix calculations 
using the fiducial model of Section IV B; future versions 
of the code will cover a wider region of parameter space. 

It is important to realize that the \ 2 statistic above 
is not necessarily the most useful or informative way of 
assessing the accuracy of the approximate curves. One 
reason is that for some spectra, the major contribution 
to the x 2 values comes from a certain range in I. For 
example, for BB tensors in the second line of Table II, 
it is only the C^s for £ > 1100 that display significant 
deviations from the exact curve. However, ignoring the 
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range 2000 < I < 3000 in our Fisher matrix calculations 
only affects our forecasts at the ~10% level, so in this 
case the sum in Eq. (1) runs over angular scales that are 
not significant for this application. 

Also, the accuracy of approximate curves is highly de- 
pendent on the spacing and extent of the pre-computed 
grid of models. Thus, in regions of the grid where the 
approximate curves are less accurate (our tests revealed 
£ < 0.15 to be one such region, if the grid spacing in £ is 
0.05 as in Table I), the results can be improved by refin- 
ing the grid and introducing more pre-computed spectra 
in these regions. Such refinements of the grid entail a 
trade-off between accuracy of the output and memory 
required at runtime, but the latter is not expected to 
present significant obstacles to obtaining usable spectra 
from the code. Future versions of the code, making use 
of larger and finer grids, will be able to produce more 
accurate curves where necessary. 

Another reason to be cautious in interpreting the ta- 
bles above is that they compare the accuracy of the ap- 
proximate curves to the variation in the "exact" curves 
from CMBACT (run at a high number of realizations), 
but it is likely untrue that Xred < 1 is necessary for ef- 
fective use of the spectra, given the shortcomings of the 
model for strings on which CMBACT is based. In addi- 
tion, even the "exact" curves will have some deviation, 
so a comparison to these curves is not completely equiva- 
lent to a comparison to the idealized limiting curves that 
are ostensibly approached as the number of realizations 
is taken to be arbitrarily large. 

Thus, it is recommended that the accuracy of String- 
Fast's output should be assessed on an application-by- 
application basis, with due consideration of the precision 
that might be desired in various ranges of £. In using 
StringFast to compute numerical derivatives of Cgs with 
respect to the parameters of the model, we assert that 
the accuracy of the approximate curves near our chosen 
fiducial model (the second line of Table II) is sufficient, 
since small and non-systematic errors in the curves will 
in general be washed out by the sum over £ in the com- 
putation of the Fisher matrix [Eq. (3)] . 



IV. METHODS 



where the Fisher matrix F for CMB observations is given 
by 



Fisher Matrices 



We obtain projected uncertainties on determinations 
of the string model parameters from the Cramer-Rao in- 
equality [51]: if Xi is an unbiased estimator for a param- 
eter A;, then 



X,Y 



(3) 



with X, Y standing for TT, EE, BB, or TE, and the co- 
variance matrix Cov incorporating signal and noise vari- 
ance on measurements of the Cgs (see, e.g., Ref. [52]). 
Thus, we take \J (F~ 1 ) ii as the expected uncertainty on 
a measurement of parameter A^. 

We calculate the derivatives in Eq. (3) numerically 
via finite-differencing, using StringFast (with the grid 
described in Table I) to calculate the string CgS and 
the public Boltzmann code CAMB [53] for the inflation- 
ary Cgs. Gravitational lensing of inflationary E-modes 
into B-modes is accomplished by CAMB; in the forecasts 
where lensing is not removed from the final signal, the 
lensed B-modcs are taken as contributing to the experi- 
mental noise in the covariance matrix. Since the E-modcs 
generated by strings arc ^2 orders of magnitude weaker 
than those generated by inflation (see the bottom panel 
of Fig. 2), lensing of these modes is ignored in our calcu- 
lation. 

We investigate the prospects for constraining the prop- 
erties of strings using the Planck satellite, and also using 
some idealized, cosmic variance-limited measurement of 
CMB polarization, which we refer to as "FuturePol" for 
brevity. Since the main constraints on strings will come 
from B-modes of polarization, this allows us to obtain the 
most optimistic estimate of how precisely we can hope to 
detect strings via their signatures in the CMB. 

For Planck, we use characteristics of the lowest six fre- 
quency channels as presented in [24] , combining the noise 
from all six channels via 



N? 



E 



■ exp 



;hi2 



(4) 



Var A, > F 



(2) 



where c runs over the channels, Z denotes either temper- 
ature or polarization, C is the beam full-width-half-max, 
and is the experimental noise per pixel (see, e.g., 
Ref. [54]). We use all four types of C e s (TT, EE, BB, 
and TE) in the Fisher matrix. 

For FuturePol, we assume the ideal case of no ex- 
perimental noise. We utilize only EE and BB, mo- 
tivated by the idea that future polarization measure- 
ments will surpass the sensitivity threshold of tempera- 
ture anisotropies, because of the lack of significant small- 
scale foregrounds for polarization signals. However, in 
order to account for the fact that an experiment like Fu- 
turePol would already have the Planck data at its dis- 
posal when the observations are analyzed, we use the 
constraints from the Planck Fisher calculation as Gaus- 
sian priors on each parameter. This is accomplished by 
adding 1/c 2 to the corresponding off-diagonal element 
of the FuturePol Fisher matrix, where a 1 is the forecast 
experimental variance on each parameter obtained from 
F- 1 for Planck [55]. 
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G/x only 
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FuturePol 
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6.40 


0.15 


0.018 


0.64 


0.073 


0.64 
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5.98 


0.35 


0.064 


12.5 


0.52 


59.8 


3.41 



TABLE III: Forecast (la) uncertainties on estimations of G/t 
from data, with no removal of lensing from the B-mode signal. 
The titles "G/x only" and "G/x + 3 others" refer to the string 
model parameters that are left free in the Fisher matrix. For 
the latter (more realistic) case, a non-zero value of G/x cannot 
be detected at better than 2a confidence if G/x < 6.4 x 10 -8 . 
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Planck 


FuturePol 


Planck 
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6.40 


0.15 


0.0038 


0.62 


0.028 


0.64 


1.14 


0.0021 


5.74 


0.0070 


0.064 


11.4 


0.018 


57.4 


0.039 



TABLE IV: As Table III, but with 100% efficiency in re- 
moving the lensing B-mode signal. Even without contami- 
nation from lensing, in the "G/x + 3 others" case G/x cannot 
be distinguished from zero at better than 2a confidence if 
G/x < 6.4 x 10" 9 . 
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FIG. 3: Prediction ellipses for string parameters measured by 
FuturePol, with G/x = 6.4 x 10~ 7 and no lens cleaning. Dark 
green (inner) regions denote forecasts at 68% C.L., while light 
green (outer) regions correspond to 95% C.L. While some de- 
generacies are present, they are not severe enough to impede 
accurate estimation of parameter values from future data. 



RESULTS 



B. Fiducial Model 

For the fiducial model we use in our forecasts (which 
then determines the point in parameter space where 
Eq. (3) is evaluated), we choose model parameters that 
were found in Ref. [34] to give spectra that agree well 
with the results of field-theoretic simulations for an AH 
string model [35]: a = 1.05, v = 0.40, and £ = 0.35. 

For the string tension G/x, we forecast for three values: 
the current observational maximum for AH strings, cal- 
culated from current CMB data and galaxy surveys [34] 
to be 6.4 x 10 -7 ; and reductions of the tension by one and 
two orders of magnitude, 6.4 x 10~ 8 and 6.4 x 10~ 9 . The 
first case allows us to investigate the "best-case scenario" 
for the detection of strings, in which the string contribu- 
tion to the CMB is as high as possible, while the other 
cases illuminate our prospects for detecting less-energetic 
strings. The background cosmology is described using six 
parameters (fl^h 2 , tt c h 2 , Hq, A^j, n s , and r) with values 
taken from the WMAP seven-year data analysis [50]. 

Previous work on strings (e.g. [32]) has not taken a, 
V, and £ as measurable quantities, instead making a pri- 
ori assumptions about the corresponding properties and 
therefore the shape of the resulting power spectrum. 
Thus, to compare the previous approach with the current 
work, we perform forecasts in two ways: first, letting a, v, 
and £ be free, adjustable parameters (with corresponding 
elements in the Fisher matrix), and second, fixing them 
at their AH values and not including them in our Fisher 
matrix calculations. 



The forecast la uncertainties on estimations of G/x 
from future data are shown in Table III, where the lens- 
ing B-mode signal is used as additional noise in the co- 
variance matrix, and Table IV, where we assume that the 
entire lensing signal can be subtracted. There are several 
notable features of these results. First, moving from the 
"G/x only" model to the (more realistic) "G/x + 3 others" 
model increases the uncertainty on estimates of G/x by a 
factor ~2-7. Since the three other string parameters can 
be estimated from CMB observations, it is therefore im- 
portant that they be allowed to vary in a thorough search 
for signatures of strings — fixing their values a priori could 
lead to erroneously precise determinations of the string 
tension. Fig. 3 shows anticipated correlations between 
measurements of the parameters; several parameters are 
weakly degenerate due to certain similar effects they have 
on CMB spectra (e.g. scaling the overall amplitude). 

Second, if the "G/x + 3 others" analysis is used, and if 
the string tension is equal to the AH observational max- 
imum, then Planck should have sufficient precision in its 
CMB measurements to determine that the string tension 
is non-zero at high statistical significance. Thus, it will 
be able to make a definitive detection of the presence 
of strings in the Universe (provided they are described 
by the AH model we forecast for here). However, if the 
string tension is just one order of magnitude below the 
observational maximum, then the efficiency with which 
the effects of lensing can be removed from the B-mode 
signal impacts the constraints on G/x. With 100% re- 
moval, FuturePol can measure the tension very precisely, 
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while with no removal, it cannot discriminate between 
the presence (G/i ^ 0) and absence (G/i = 0) of strings 
at better than -~2cr confidence. For G/i two orders of 
magnitude below its upper bound, even FuturePol with 
100% lensing removal cannot do any better than a ^2er 
detection of a non-zero G/i through its measurement of 
the G^s. 

We can compare this result to a recent similar 
study [32], which asserts that AH-generated strings can 
be detected by CMBPol [56] at 3cr confidence if G/i ~ 
9 x 10~ 8 . That work uses the cosmic variance of the weak 
lensing B-modcs as additional Gaussian noise, while we 
use either the entire lensing signal as noise or assume that 
lensing can be completely removed in the final analysis. 
We agree that after assuming specific shapes for the Ct 
curves, so that the only freedom left is to scale the over- 
all amplitudes via G/i, CMBPol or our FuturePol will be 
able to detect a nonzero string tension at high confidence, 
even without removal of any contamination from lensing 
(see the "G/i only" section of Table III). 

However, we have shown that there are two factors 
that can strongly impact the precision with which we can 
measure G/x: the efficiency of lensing removal, and the 
way the properties of strings are allowed to change the 
shape of the CMB power spectrum. We reiterate that it is 
preferable to relinquish the assumption of a specific shape 
for the spectrum (an effective fixing of a, v, and £) and 
allow other properties of strings to be measured from ob- 
servations, in which case strings with G/i < 7 x 10~ 9 can- 
not be reliably constrained even by a noiseless, lcnsing- 
free measurement of polarization power spectra. Thus, 
theories that predict a wide possible range of orders of 
magnitude for G/i will require observational probes other 
than the CMB Cis to be verified or ruled out completely. 

VI. CONCLUSIONS 

Our primary goal in this paper has been to predict how 
accurately the tension of cosmic strings can be measured 
by upcoming measurements of CMB polarization (and 
temperature, to a lesser extent). To do this, we have de- 
veloped a new scheme to quickly calculate the Cgs from a 
network of strings with certain specified properties (ten- 
sion, mean velocity, correlation length, and "wiggliness" ) , 
and implemented this scheme in a modular code that we 
call StringFast. Our calculations have the advantage of 



incorporating these properties as adjustable parameters 
of the model for strings, as opposed to fixing these prop- 
erties a priori by assuming a specific shape for the power 
spectrum generated by strings, which is often done in 
the literature. Using a Fisher matrix approach, we have 
estimated how confidently we can constrain G/i with up- 
coming CMB data, and also how much this constraint is 
weakened by adding additional physically-motivated free 
parameters to the numerical description of strings that 
is employed. 

Our results highlight the importance of prior assump- 
tions in the search for cosmic strings in observational 
data. It is advisable to allow for more general string 
properties than those of any one field-theoretic model 
or numerical implementation, and thereby incorporate 
the uncertainty inherent in computational approaches to 
modelling strings — but, doing so also weakens our ability 
to detect strings via their effects on the CMB. Our fore- 
casts also indicate that theories that generate strings with 
G/i <C 10~ 9 do not have strong prospects for verification 
via CMB power spectra, and thus we may need to look to 
gravitational wave signals or other signatures for empiri- 
cal tests. Note, however, that our analysis only applies to 
field-theoretic strings and not to large-scale supcrstrings; 
see Ref. [37] for a similar recent study of superstring net- 
works with multiple tensions and Y-junctions. 

With the impending release of data from Planck, and 
with other upcoming probes of polarization, we will soon 
be able to place new, tighter constraints on the energy 
scale and properties of cosmic strings. While the search 
for strings is not the primary aim of these experiments, 
their data can and should be used for an attempted detec- 
tion (or null detection) of strings, either of which would 
provide important insights into the physics of the early 
Universe. 
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